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Abstract 

The results of a long-established investigation into pupal transpiration are used as a rudi- 
mentary data set. These data are then generalised to all temperatures and humidities by 
invoking the property of multiplicative separability, as well as by converting established 
relationships in terms of constant humidity at fixed temperature, to alternatives in terms of 
a calculated, water loss rate. In this way a formulation which is a series of very simple, 
first order, ordinary differential equations is devised. The model is extended to include 
a variety of Glossina species using their relative surface areas, their relative pupal and 
puparial loss rates and their different 4th instar excretions. The resulting computational 
model calculates total, pupal water loss, consequent mortality and emergence. Remaining 
fat reserves are a more tenuous result. 

The model suggests that, while conventional wisdom is often correct in dismissing vari- 
ability in transpiration-related pupal mortality as insignificant, the effects of transpiration 
can be profound under adverse conditions and for some species, in general. The model 
demonstrates how two gender effects, the more significant one at the drier extremes of 
tsetse fiy habitat, might arise. 



Keywords: pupal water loss; transpiration; dehydration; mortality; emergence; tsetse; Glossina 
morsitans; Westwood; trypanosomiasis; gambiense; rhodesiense; sleeping sickness; diptera; 
glossinidae. 



1 Introduction 



Early stage mortality is considered to be, by far, the most significant effect in any model of 
tsetse population dynamics (HARGROVE (Si) and pupal dehydration appears to represent the 
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most challenging aspect of modelling it. The implications of pupal dehydration are far greater 
than pupal emergence and mortality alone. Water loss continues after eclosion up until the 
moment the teneral has its first meal. To give some idea of relative importance, it can be 
argued that, while teneral water loss rates are generally 20 times puparial rates and 100 times 
pupal rates (comparing Bursell [1] with BURSELL |l2l data), puparial rates generally prevail 
6 times longer and pupal rates 24 times longer than teneral rates. Thus, any teneral that dies 
of dehydration could be said to be, at very least, 35% as likely to have died as a result of 
pupal water loss. If it is further specified that the teneral was of average age, the figure is 
closer to 51%. Water loss during the pupal phase can decide the fate of the teneral. Combined 
dehydration and fat loss are thought to culminate in massive teneral mortality. Pupal and teneral 
mortality rates are crucial in deciding the viability of any Tsetse population. The ultimate effect 
of cumulative water loss on a given cohort is therefore likely to be best assessed in terms of 
the proportion which have sufficient reserves to achieve their first feed. The vastly different 
dynamics of water loss during the pupal and teneral phases, however, afford pupal water loss 
the status of a topic in its own right. 

This work is almost exclusively based on the findings of one experimentalist with all the haz- 
ards implied. In 1958, the late E. Bursell published the results of his experiments on pupal 
water loss. Today, in an age of mainstream computing, that work turns out to be a somewhat 
tantalizing, scientific riddle. As one might imagine, Bursell's results are not of much use in the 
form in which they were presented. Most of the work was carried out for steady humidities at 
24.7°C. The main challenge to exploiting it for the purposes of a computational model, lies in 
generalising the results to all temperatures and humidities. That challenge could be said to be 
in three, very specific respects: A function for transpiration, the historical conditioning of sensu 
strictu pupal transpiration, then linking water loss to observed pupal emergence. Two further 
obstacles arise in the form of, firstly, extending the Glossina morsitans-hased model to the rest 
of the Glossina genus and, secondly, resolving the dependence of emergence on humidity. 

The transpiration data were assigned to one of three categories for the purposes of this work: 
Temperature dependent data, humidity dependent data and time dependent data. One observes 
a certain amount of corroboration between points on the respective curves. Some of this cor- 
roboration is demanded, for example, at intersections, however, other of it comes as a pleasant 
surprise. The time dependent data is a case in point. Time dependent transpiration would 
appear to be nothing more than an alternation between pupal and puparial rates. The final for- 
mulation, hence solution to the problem, is predicated on five major assumptions. In addition 
to the assumptions explicitly stated and explored, two others are taken for granted. The first 
is that it is assumed that Bursell yj is comprehensive, to the extent that it encapsulates all 
salient aspects of pupal water loss. The second is that there is no transpirational water loss at 
dewpoint. The problem is then reduced to a series of first order, ordinary differential equations 
for water loss. 

Although these equations are extremely simple, they are both numerous and voluminous. 
The relevant domains of applicability are also non-trivial. This renders preferred integra- 
tion schemes, such as the fourth order accurate Runge-Kutta-Fehlberg (R-K-F-4-5) method 
slightly impractical. Since the problem is not intractably large, expedience takes precedence 
over taste and the more pedestrian Euler's method is the preferred integration technique. The re- 
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suiting problem becomes a computational one, rather than a mathematical one. A least squares 
fit, Newton's method and a half-interval search are the only other techniques employed, the 
mathematics in the model being nothing more than utility. 

This model is an experimental model. "Experimental" in the sense that is based on data gleened 
mostly from pre-1960s-published graphs and it relies almost exclusively on the work of one 
experimentalist. It assumes that issues such as inferior quality pupae, differing puparial dura- 
tions and the shortage of statistically significant data can be rectified at a later stage. Extending 
a G. Morsitans-hased model to other species is work that can best be described as exploratory. 
The results at the end make it an interesting and justifiable exercise, nonetheless. 

2 Generalising Scant Transpiration Data to a Function of 
Variable Humidity and Temperature 

BURSELL |[T1 obtained one set of data points for variable temperature (at 0% r.h.), one for 
variable humidity (at 24.7 ± 3°C) and another for temporal dependence (at 0% r.h. and 24.7 
± 3°C), during his investigations into pupal water loss. Yet a fourth set of data points can be 
inferred by reason. One expects no transpirational water loss at dewpoint, regardless of the 
temperature. Although of some assistance, the challenge, nonetheless, remains: How does one 
generalise these data to all temperatures and humidities? Fortunately, enough of the aforemen- 
tioned data exist to suggest that any transpiration function is not only continuous, it is also 
surprisingly simple and smooth; monotonic, in fact. 

Assumption 1 Transpiration rate is a multiplicatively, separable function of humidity and 
temperature. Put succinctly, if ^ is the transpiration rate, then there exist two functions <j) 
and 9, dependent exclusively on humidity and temperature respectively, so that 



in which h denotes humidity and T, the temperature. 

Just how reasonable is this assumption? Certainly it is consistent with, and replicates, the 
fourth, inferred set of data points entertained above. The perceived wisdom is that the region 
of interest is between 16°C and 32''C (due, in part, to other causes of mortality). For the "H" of 
data which exists across the humidity-temperature domain one reasonably expects rates to be 
bounded by the wet and dry end data, furthermore, to be close to monotonic. One would also 
expect any unusual, capricious behaviour, or even failure in the waterproofing, to manifest itself 
in dry air. The dry-air, temperature-dependent data set is, fortunately, reasonably complete 
and suggestive of behaviour which is smooth, monotonic and simple (pure exponential, in this 
case). Thus, in the very likely event that water loss rates are not multiplicatively seperable, 
multiplicative seperability should not be a bad substitute. 



dk 
'dt 



{h,T) 



(1) 
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Of course, one can never be sure in these matters. Anything is possible. As much as someone 
who models with data in the ideal format of a grid is ultimately ignorant of the behaviour 
between grid points, one is faced here with the same possibility of some magic combination of 
humidity and temperature. One simply doesn't know; one can only surmise. Engineers make 
the same assumption under what are sometimes, seemingly, a lot less favourable circumstances. 
Their justification? It works. 

The multiplicative separability assumption fulfills a secondary purpose in that it facilitates the 
following legerdemain 

_ 'iiw,h,2A.7) e{T) 
dt^ ' ' ' ~ ^(24.7) 

the need for which will become clear further on (w is a total historical water loss). 



2.1 The Dependence of Transpiration on Temperature 

Two types of transpirational temperature dependence are recognised in keeping with the mul- 
tiplicative separability assumption. The first is a puparial type and the second is a pupal type. 
Both are exponential in nature. 



2.1.1 For the Puparium 



Certainly so far as transpiration rates are concerned, for the puparium it is a case of the worst 
first. Puparial transpiration rates oscillate wildly over time and a fairly substantial difference 
in data which should corroborate is documented. An average of the relevant data points from 
BURSELL Q Figures 2a, 2b, 3 and 8a (those at 0% r.h. and 24.7°C) were accordingly used 
to adjust an exponential fit to the Figure 8a data upward. The basis for this decision was 
threefold. Firstly, the Figure 8a data were read from a logarithmic curve with a consequently 
greater, implied possibility of error. Secondly, the humidity data were more comprehensively 
presented. Thirdly, the study was overridingly an in-depth study of the effects of humidity, 
strongly suggestive of an overall greater attention to detail and accuracy pertaining to humidity 
dependence. The following function was the result, 

Wium(r) = e°-^^02^«^-9-^22°^ + 0.000354783, (2) 

the units of which are G. morsitans, initial pupal masses per hour. The assymptotic standard 
errors' in fitting the constants for the exponential power were 6.8% and 3.32% respectively. 
The sum of squares of residuals' was 2.98448x10"^. 



It is for consistency with this information that no attempt to guess the number of significant figures has been 
made. 
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Figure 1 : The fit to the Bursell Figure 8a data were adjusted upward. 



Of course, one might just as well have used the aforementioned graphs to argue the case for an 
upward adjustment of 0.00039619. There are many alternatives, however, a decision had to be 
made and so 'modeller's licence' was invoked to make the choice which best takes cognizance 
of the, as yet unused, Bursell Q Figure 1 data. 



2.1.2 For the Pupa 



Relevant data points from Bursell ^ Figures 5b and 6 (those at 0% r.h. and 24.7°C) were 
added to the Figure 8b data of the same author. The following fit was obtained, 

^pupa(T) = e'-''''''^-''-'''\ (3) 

the units of which are G. morsitans, initial pupal masses per hour. The assymptotic standard 
errors' in fitting the above constants were 4.597% and 2.651% respectively. The sum of squares 
of residuals' was 3.62324x10"^. 



Note that, strictly speaking, the BURSELL fT\ Figure 8a data only pertains to the first day of 
the puparial duration. This will be of relevance in devising a temporal dependence. 
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Figure 2: The fit to the Bursell Figures 8b, 5 and 6 data. 



2.2 The Dependence of Transpiration on Humidity 

Two basic types of transpirational humidity dependence are recognised in keeping with the 
multiplicative separability assumption. The first is a puparial type and the second is a pupal 
type. Both are dimensionless. 



2.2.1 For the Puparium 

A cursory inspection of BuRSELL Hj Figures 2a and 2b leads to the deduction that 

, _ 100 -h 

The dependence of transpiration on humidity is linear. 
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2.2.2 For the Pupa 

Pupal transpiration is not as straightforward. While the above relation may prevail during the 
initial transition down to pupal rates, one dependent on both total historical water loss and 
humidity is ultimately required. Historical conditioning, what one might term 'drought hard- 
ening', alternatively depletion, is a phenomenon which pertains to Tsetse pupae. During the 
sensu strictu pupal phase, transpiration becomes conditioned by the temperature and humidity 
which prevailed during the early stages (3rd and 4th instars inclusive). Present transpiration is 
conditioned by the recent past, in addition to the prevailing humidity and temperature. 

Relevant data is that published in Figures 5a, 5b and 6 of BURSELL While it was, no doubt, 
acquired with a different purpose to the present one in mind and certainly proves a point, it is, 
of little use as presented. The problem is that historical humidities are steady, furthermore, the 
data were obtained at 24.7"C. 



Assumption 2 The transpiration rate, conditioned by a given historical water loss, is the 
same as the transpiration rate conditioned by an historically steady humidity, at 24.7^ C, 
which produced an equivalent total water loss. 



In other words, a conversion of the independent variable, historically-steady-humidity-at- 
24.7°C, to an associated total water loss is inferred. Transpiration can then be re-expressed as 
a function of historical water loss. An historical conditioning of the pupa which is dependent 
on historically-non-steady variables is devised in this way. How reasonable is the assumption? 
Do different histories in temperature and humidity, which produce the same water loss, imply 
the same historical conditioning? If not, there are additional historical effects that have never 
been detected. 

The BuRSELL 01 Figures 5a, 5b and 6 data can be interpretted as transections through a surface 
which intersect at their ends, once their dependence on humidity has been converted to one of 
total historical water loss. They suggest a very simple surface, one which appears to be of no 
higher order than bi-quadratic, by inspection. It was therefore decided to fit a bi-quadratic 
surface to the historically conditioned transpiration data, that is 



using the method of least squares, 'on the fly' so-to-speak. The results were pleasing in that the 
surfaces retained their fundamental character, even for fictitious, negatively-large humidities, 
the importance of which will become apparent further on. 



dk 
'dt 



pupa 



{w,h,24:.7) 



ci + C2h + C3W + Ciwh + c^h^ + ceui^, 
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Figure 3: The hourly transpiration rate as a function of humidity and water lost during the first 
8/30 of the puparial duration (for G. morsitans). 



The surface is, nonetheless, a surface of transpiration when, instead, the humidity dependence 
is sought. The following legerdemain 



dk 
dt 



>,/i,24.7) 



is based on the multiplicative separability assumption. Thus, the humidity dependence 

/ , - Ci + C2h + C3W + C4Wh + Cs/i^ + CgW^ 

is obtained. 



3 The Dependence of Transpiration on Time 

Time dependent transpiration would appear to be nothing more than an alternation between 
pupal and puparial rates. Only the following stages of the puparial duration, r, were deemed 
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worthy of any time-dependent modelling, the vagaries of which were entirely reduced to visu- 
ally fitting four straight lines to the as-yet-unused BURSELL |[T| Figure 1 data. 



The Period to 



During this period there is an adjustment from the puparial rate down to the pupal rate so that 
the total transpiration rate can be approximated as 

dk _ dk _ 27 jt-^r) / dk _ ^ ^ 

dt dt puparium 31 — V puparium dt pupa J 

in which t is time, 

dk dk' 

J , 0puparium (^) ^puparium (-^) and , 0puparium(^)^pupa(-^) ; 

(It puparium Clt pupa 

^pupa t'eing a temporary, or transitional transpiration rate; one defined shortly prior to that for 
which the dependence on historical water loss is known. 

The Period joTto-^r 

This period represents final adjustment down to the pupal rate. The equation 

dk _ dk _27 fdk _ ^ 3 - ^t) /dk _ dk' 

dt dt puparium 31 y dt pupaiium dt pupa J 31^^T — i^^) V puparium dt pupa 

(8) 

was used. 



The Period to |§r 

Transpiration during this phase is predominantly at the pupal rate. A small component of loss 
at puparial-rates increases linearly with time. The resulting combination was deemed to be 

dk _ dk ^1 {t-i^r) /dk dk \ 

dt dt pupa 30 — iro^) \ '^^ P"^?^™^ dt pupa J 

in which 

dk ^ f (^;,/i,24.7)gpupa(r) 

(it pupa ^pupa(24.7) 

based on equation [51 
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The Period |§r to fgr 

Transpiration begins its return to puparial rates during the pharate adult phase and is modelled 
by 

dk _ dk ^}_fdk _dk \ 7 (t-fr) /dk _dk \ 

dt dt pupa 30 y C?t pupai'ium dt pupa J 30 (^''^ — l^^) \ puparium dt pupa J 

(10) 

in this work. 



4 Extending the Model to Other Species 

It is generally suspected that the Glossina genus derives from a common, tropical, rain-forest 
dwelling ancestor, adjusted to moist, warm climates. One might therefore also suspect that 
all tsetse species actively pursue a strategy to minimise water loss for the majority of modem 
habitats and have hydrational mechanisms preventative of desiccation. It is generally accepted 
that most of the genus is not well adapted to arid environments GLASGOW [i4|. The challenge 
to pupae, indeed the major threat, is dehydration. 

Assumption 3 The hydrational mechanisms and water management strategies of the ma- 
jority of Tsetse fly species differ only with respect to relative pupal surface area, relative 
puparial loss rates, relative pupal loss rates, the different amounts excreted during the 4th 
instar and initial reserves. 

Water loss rates for the puparium and pupa, Ppuparium and Ppupa, respectively, have been measured 
for a number of species and are tabulated in Bursell [[T]. Permeability is dependent on 
pressure and is quoted in units of mg h^^cm~^(mm Hg)~^. No variation with pronounced 
varaiation in temperature and humidity is indicated and it is of some comfort that the conversion 
of the G. morsitans model to other species involves relative rates. 

Surface area data is likewise available. The same surface area is used for both the puparium 
and pupa in this work, the justification being that the puparial exuviae render the puparium 
marginally bigger while the pupal surface is not as regular. 

4.1 For the Puparium 

A dimensionless, species conversion factor for puparial transpiration rates can be defined as 
follows 

e Ppuparium •5 species 

'^puparium ^ i 

Pmorsitans puparium "^morsitans 
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in which Ppuparium is the rate of water loss for the species in question, Pmorsitans puparium is the 
equivalent water loss for G. morsitans, Sspedes is the puparial surface area for the species in 
question and s^orsitans is the equivalent, G. morsitans surface area. This factor enables the 
puparial transpiration rate for another species to be calculated from G. morsitans values. Note 
that the unit is still in G. morsitans initial pupal masses (31mg) per hour. Actual values of 
^puparium fof tcu different species are tabulated in Table [T] 



Group 


Species 


'^puparium 


'^pupa 
l^iur lllllllllla^ 


^pupa 
l^iUr lllaAlllla^ 


morsitans 


austeni 


1.60 


0.712 


0.723 




morsitans 


1 


1 


1 




pallidipes 


1.50 


1.24 


1.31 




submorsitans 


2.44 


0.950 






swynnertoni 


0.830 


0.869 


0.892 


palpalis 


palpalis 


2.54 


1.41 


1.36 




tachinoides 


1.26 


0.743 




fusca 


brevipalpis 


10.3 


4.57 


3.06 




fuscipleuris 


8.84 


4.45 


3.16 




longipennis 


3.62 


2.45 


2.30 



Table 1: Species conversion factors for the model calculated from data, ultimately sourced 
from Buxton and Lewis [SJ, presented in Bursell [jy. 



4.2 For the Pupa 

A dimensionless, species conversion factor for pupal transpiration rates can be defined as fol- 
lows 

r Ppupa '5 species 

'^pupa ^ 

Pmorsitans pupa "^morsitans 
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in which ppupa is the rate of water loss for the species in question, Pmorsitans pupa is the equivalent 
water loss for G. morsitans, Sgpedes is the pupal surface area for the species in question and 
Smorsitans is the equivalent, G. morsitans surface area. This factor enables the puparial transpira- 
tion rate for another species to be calculated from G. morsitans values. Note that the unit is still 
in G. morsitans initial pupal masses (31mg) per hour. Actual values of 5pupa for ten different 
species are tabulated in Table [B 

On the face of it. Assumption [3] is certainly the most tenuous. How valid is it? Does such a 
simplistic approach work? Very little data is available for other species, however, 3rd instar, 
puparial loss rates for both Glossina brevipalpis and Glossina palpalis are known. Conversion 
of 3rd instar, G. Morsitans-model, transpiration rate values to G. brevipalpis and G. palpalis 
values, yielded errors of 6% and 10% respectively. The Spupn for pupal maxima and minima 
in Table \T\ are, furthermore, remarkably similar (for all except G. brevipalpis and Glossina 
fuscipleuris). This is very encouraging and suggestive of a similar slope in the transpirational 
time dependence for the various species. The suggestion for G. brevipalpis and G. fuscipleuris, 
however, is that As sumption [3] could possibly be captious. 

What are the implications? It means that knowing only the appropriate G. morsitans (j) and 
9 is adequate. The only questions pertaining to species conversion which remain are whether 
the temporal interplay between pupal and puparial transpiration rates is the same for the entire 
puparial duration and whether the historical conditioning is the same; not withstanding some 
difference in strategy a la the difference in 4th instar excretions ^ 



5 The Resulting Model for Pupal Water Loss 

Taking into account one, further formula (that for puparial duration), results in a model. 



5.1 The Puparial Duration 

The puparial duration in days, r, is calculated according to the formula 



in which k = 0.057 ± 0.001, a = 5.5 ± 0.2 and b = -0.25 ± 0.01, for females and n = 
0.053 ± 0.001, a = 5.3 ± 0.2 and b = -0.24 ± 0.01, for males (Phelps and Burrows BH 
modified by HARGROVE [|6l). The puparial durations of all species, with the exception of G. 
brevipalpis, are thought to lie within 10% of the value predicted by this formula Parker [|8l. 
G. brevipalpis takes a little longer. 

Newton's method is used to solve for a puparial duration based on the daily average, which is, 
of course, dependent on itself. The same applies for the various fractions of puparial duration. 



Although, in that case, the associated water loss is too small to be of any real consequence. 
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For this reason the notation ts. is adopted, where 

floors T r \ 

\ tot/ 

{ts. - floor {t^ }) ^ r(Tj^^ ^^^^|^^ |^^) + ^ r{T,,y ^) 

i=l 



'^30 J- ^ 

30 



it being the average x t over the time interval (0, rj_], for a specified r. 



5.2 The Governing Equations 

Collecting together all prior observations and thoughts gives rise to the following series of first 
order, ordinary differential equations. Note that what would otherwise have been a unit of 
31mg h^^ can be replaced by a dimensionless 

^morsitans • •.• i i 24Ti 

initial pupal masses x 



mspecies 1 puparial duration 

for the usual reasons, where m,T,orsitans and mspecies are the initial pupal masses of G. morsitans 
(31mg) and the species in question respectively. In other words, a dimensionless rate unit of 
initial-pupal-masses per puparial-duration is preferred. 



The Period to 



30 



The water loss rate for the greater part of the third and fourth instars is at the puparial rate 
(obtained by the substitution of equation [2] and equation |4] into equation [1]). Generalising the 
resulting expression to all species and writing the equation in dimensionless form results in 

^ _ |- 0.110268T-9.92201 , Q 000354703') ~ ^ Ppupaiium ^species ,^2\ 

xww /^morsitans pupanum ''morsitans 



The Period r j_ to t_6_ 

30 30 

During this period there is an adjustment from the puparial rate down to the pupal rate dictated 
by equation|71 Generalising equation|7]to all species and writing the equation in dimensionless 
form results in the expression 

dk 
H 



^g0.110268T-9.92201 ^ Q00354783) I 1 _ ^ j Pv^V^nn^ 

31[-j-g T4I / Pmorsitans puparium 



+ e' 



0.161691T~12.9591 2''' ^ 30 ) Ppupa 



31|-|-g T4) Prnorsitans pupa 

V 30 30 / 



100 h Sspecies (13) 

100 

-^morsitans 



in which t' is a developmental 'time', t = t_l. 
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The Period Te_ to ts_ 

30 30 



During this period there is a final adjustment from the puparial rate down to the pupal rate 
dictated by equation [8l Generalising equation [8] to all species and writing the equation in 
dimensionless form results in 



dk 
'dt 



.0.110268T-9.92201 ^ 0.000354783) I - - - 

^ ' 31 31 



+ e 



,0.161691T-12.9591 



27 3 
31 ^ 31 



t' - TJL 

30 



T8_ — TG_ 

30 30 




Ppuparium 



Pmorsitans puparium 
100 hy <5species 



100 s„ 



(14) 



Excretion 

If water loss is sufficiently low during the first of the puparial duration, cognizance must be 
taken of the small amount excreted. In this unlikely scenario the formula 

h 

7 , ' ^3rd instar / \ 

k = Xo -\ [Xl — X2) 

100 ^ ' 

was implemented, for want of any better wisdom. The total water loss during the 3rd and 4th 
instars, in the event of dewpoint prevailing for the former, is x^. In the event of 0% relative 
humidity prevailing for the 3rd instar, the amount is X2. The only 4th instar excretion data 
known to exist is that for G. morsitans, G. palpalis and G. brevipalpis. This lack of information 
is a minor obstacle as the excretions are generally small and only relevant for humidities close 
to dewpoint. 



The Period t_8_ to 

30 30 



Transpiration during this phase is predominantly at the pupal rate. There is also deemed to be 
a small component of loss at puparial-rates, which increases linearly with time and which is 
included for good measure. Generalising equation |9] to all species and writing the equation in 
dimensionless form results in 



dk 
'dt 



g0.161691(T-24.7)(^^ + C^h + C^W + C,wh + C^h^ + C^W^) 1 



t' - 

30 



30 



r 25 

30 



30 



^0.110268T-9. 92201 



Pmorsitans pupa 
1 



0.000354783) 



100 - h 
100 



t' - Ta_ 

30 



Ppupaiium 



30 



T25 — T8_ 

30 30 



Pmorsitans puparium 



-"species 



(15) 
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The Period T25 to 

30 30 



Transpiration begins its return to puparial rates during the pharate adult phase. Generalising 
equation [10] to all species and writing the equation in dimensionless form results in 



dk 
Itt 



30.161691{T-24.7) 



Ppupa 
Pmorsitans pupa 



,17 v^'-^i, 

X + 



30 30 



Ci + 02/1 + C3W + C4,wh + Cs/i^ + CqIV^) I — 



,o.iio268r-9.922oi ^ 0.000354783) 



1 29 


7 


(''- 


T2b ] 
30 / 


1 30 


30 




- T25 ) 






\ 30 


30 / 


-h 









30 30 , 



Ppuparium 



Pmorsitans puparium 



100 



^species 



(16) 



The Period to ri 



There is a return to puparial rates shortly before eclosion and equation [121 once again applies 



5.3 Solving the Equations 



The above rate formulae constitute a series of first order, ordinary differential equations. One 
expects the resulting function to be Lipshitz continuous over each of the developmental sub- 
stages identified, likely even a contraction. While a fourth-order-accurate Runge-Kutta- 
Fehlberg method (R-K-F-4-5) would normally be the preferred method of integration, the 
series of equations is voluminous and the relevant domains of applicability are also non-trivial. 

Euler's method is usually considered distasteful from the point of view of its error. The local 
error per step, of length At, is O(At^). Since the required number of steps is proportional to 
the global error is 0(At). This is indeed primitive. The method is, nonetheless, considered 
robust for the type of first order, ordinary differential equation to be solved. The real strength 
of Euler's method lies in its robustness at discontinuities and points of non-differentiability. 
The maximum, additional error introduced at such points is of the same order as the method's 
global error (this is easy to see). The same cannot be said for the higher order methods. The use 
of one or other of the higher order methods is still not precluded in the problem at hand, since 
the discontinuities and points at which differentiability breaks down are predictable. Using 
Euler's method, however, one has one problem to solve, whereas using one of the higher order 
methods entails solving six, seperate problems; each confined to its own respective domain of 
Lipshitz continuity, scaling etc.. 

The handicap of a poor error is easily overcome computationally. That is by using a small 
step length e.g. jj^th of a puparial duration. Two significant figures are all that is sought. 
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Since the problem is not intractably large, expedience takes precedence over taste and the more 
pedestrian Euler's method is considered the appropriate choice. 



6 Pupal Emergence and Mortality 

Two challenges arise when it comes to pupal emergence: The first is to establish some kind of 
credible relationship between the numbers of emergent and humidity. The second is, conse- 
quently, how to relate emergence to total water loss. 



6.1 What is the Relationship Between Pupal Emergence and Humidity? 

What does one make of the very rudimentary data in Figure |4l? There are insufficient data 
points for each species to perform any kind of rigorous hypothesis testing. 




20 40 60 80 100 

relative humidity / % 



Figure 4: Emergence data as presented by BuRSELL |[T1. All are at 24°C, except Tachinoides 
(30°C). 
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Since the number of pupae contributing to each data point is so low, neither the law of large 
numbers can be invoked, nor therefore, the central limit theorem applied. This does, nonethe- 
less, not necessarily preclude the use of the aforementioned in argument. Large numbers of 
pupae do exist. When doctrinaire methods fail a little thought can still go a long way. 

At the simplest level, one would expect each species to be adapted to some ideal humidity 
for which emergence is optimal. One also expects the individual pupae of each species to 
exhibit a certain amount of variation about the mean so far as size, reserves, competency of the 
integuments and so on, is concerned. Some pupae will be slightly bigger, have slightly bigger 
reserves and more competent integuments. Yet others will be slightly smaller, have slightly 
smaller reserves and less competent integuments. To be succinct, one expects emergence to be 
Gaussian. 



Assumption 4 The relationship between pupal emergence and humidity is a Gaussian 
curve, or a part thereof. 



At one extreme, one has environmentally highly specialised species with low hydrational inertia 
(e.g. Glossina austeni^ and G. brevipalpis^) , for which one expects variation over a small range 
in conditions to provide adequate data to fit the Gaussian curve. 

At the other extreme, species which exhibit massive hydrational inertia exist, such as G. longipen- 
nis and G. swynnertoni. They provide little, or no clue as to the underlying Gaussian relation- 
ship between emergence and humidity. The range of conditions, the domain, is not obviously 
suggestive of an underlying Gaussian emergence curve. All one sees is a very small, con- 
sequently flat-in-appearance, sample of the top of the curve. These species ought to have 
been investigated in terms of water loss rather than humidity. The curves for other species lie 
between these extremes. 

In the wild, there is a compounding factor in that, not only is emergence based on variation 
within a given species, it is also based on variation within the environment. A whole range of 
microclimates exist within breeding sites, some of which are compost, rot holes and soil, to 
name only a few. Seasonal variation is a further compounding factor. 

While the focus of this work is desiccation it is of interest to note that emergence also declines 
at very high humidities. As to whether drowning or some fungus is the desiderate explanation, 
it can only be speculated. 

A chi squared test is not expected to elucidate any more than visual inspection. Despite the im- 
possibility of any rigorous hypothesis testing the author maintains that the assymptotic standard 
errors obtained in Table |2] make a compelling argument for As sumption HI 



' G. austeni has low hydrational inertia by virtue of its small size (which implies a high surface area to volume 
ratio), G. brevipalpis has low hydrational inertia due to inferior waterproofing. 
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Figure 5: Species with low hydrational inertia (e.g. G. austeni and G. brevipalpis) are good in- 
dicators of the underlying Gaussian relationship between emergence and humidity. (Particular 
attention is drawn to the data for the Brevipalpis fit in Table |2] on page [191) 




20 - 



G. swynnertoni 

" " G. longipennis 

20 40 60 80 100 

relative humidity / % 

Figure 6: Species with high hydrational inertia (e.g. Glossina longipennis and Glossina swyn- 
nertoni) provide little, or no clue as to the underlying Gaussian relationship between emergence 
and humidity. 
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E{h) = ae 
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(h-73. 1591)2 


3.855 


2.077 6.401 


41.6495 


morsitans 


austeni 


lUi.OOOe 2x30.6468-i 




morsitans 


94.4792e 2x77.3495^ 


5.177 


18.45 29.07 


305.227 




pallidipes 


{h-71. 5636)2 


2.833 


6.205 9.518 


48.7588 




86.6257e 2x54.97132 

{?i-81. 1895)2 




submorsitans 


94.5092e 2x75.44742 

{h-62. 4064)2 


6.524 


37.45 57.03 


80.4127 




swynnertoni 


94.0194e 2x75.23392 

(h-78. 8419)2 


4.343 


8.65 17.19 


21.3715 


palpalis 


palpalis 


95.8732e 2x23.4835^ 


13.89 


5.412 23.94 


725.514 




tachinoides 


{?i-79. 6877)2 


11.02 


17.66 29.46 


427.406 




98.8383e 2x4o.86i62 

(Ji-84. 0199)2 


fusca 


brevipalpis 


94.0057e 2x13.6433^ 


0.5123 


0.1352 0.907 


0.268286 



Table 2: The percentage emergence, the assymptotic standard errors and the sum of squares of 
the residuals for the fit in each species. 
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Figure 7: Percentage emergence modelled as a Gaussian curve for a variety of species. All are 
at 24°C, except Tachinoides (30°C). 



The Issue of Sub-Standard, Laboratory Pupae 



In BURSELL [[II it is somewhat heuristically argued that the laboratory pupae in question were 
too small and that all emergence curves should therefore be displaced 10% to the left (the right 
in those graphs). An alternative argument based on puparial transpiration and in which the 
pupa is approximated as a spheroid, is preferred (on page|30]of the addendum). It entertains 
replacing the emergence function, E{K), with 

E { (.7) 

\ 3/100 / 
\ V 90 / 

as an alternative. 

'G. longipennis is the single exception (a straight line had to be fitted to the only two data points). 
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6.2 Survival to Emergence for Historically-Non-Constant Humidities and 
Temperatures 

Problems similar to those encountered for historical water loss, compound matters when it 
comes to the survival to eclosion for each species. Historical humidities are steady, further- 
more, the data were obtained at a constant 24"C. 



Assumption 5 The pupal emergence for a given water loss, is the same as the pupal emer- 
gence for a steady humidity at 24^ C, that produced an equivalent total water loss. 



In other words, it is assumed that emergence can be re-expressed in terms of total water loss. 
Do different histories in temperature and humidity, which produce the same water loss, imply 
the same pupal emergence, or is the amount of water present at some particular stage more 
relevant to the pupa's survival to full term? The simple answer is to refer the reader to the title, 
although this does somewhat avoid the question. 

In practice, it is far easier to convert a total computed water loss to a corresponding steady- 
humidity-at-24°C\ instead of the other way around . Either way entertains the possibility of 
fictitious, or negative, humidities. (There are always those who are apt to find this sort of thing 
vaguely disturbing, however, it should be pointed out that E{h) is just a mathematical function 
and /i, a variable. One really needs to think 'outside the box' in these matters. What was 
humidity is now not so much an artefact, just something a little more abstract.) 

The results of the water loss algorithm for 24°C and any, given set of steady humidities obvi- 
ously constitute a monotonic decline. The problem, however, is that the complete algorithm is 
relatively involved and voluminous. Under these circumstances, practical considerations and 
not rates of convergence dictate implementation of an 1/2 interval search. (The rate of conver- 
gence is not bad in this instance.) 



7 Remaining Fat Reserves 

Remaining fat reserves are a more tenuous result. Water content stays constant after the 3rd 
and 4th instars GLASGOW [4J. Although the oxidation of fat for the specific purpose of water 
production is suspected, it could not be proven at a statistically significant level (BURSELL |[ll). 
The conversion factor (by mass) is given as 

fat oxidised = 1.12 x water. 



'G. tachinoides data the exception, being at 30°C 
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8 Testing the Model 

Testing the model presents something of a challenge. The emergence data at 24°C, to a certain 
extent, provides a test for self-consistency. Transecting the following surfaces of emergence at 
24°C should come close to replicating the Gaussian curves for emergence, on page[2Ql "Close", 
since the data were adjusted as a consequence of the inferior pupae issue. 

Testing for corroboration with observed mortalities in the field is somewhat heuristic. All one 
can say is that predicted pupal mortalities due to water loss should, logically, never exceed 
any pupal mortalities observed in the field for similar conditions of humidity and temperature. 
For example, the computed G. morsitans emergence due to water loss is not lower than the 
87% obtained by HARGROVE AND WILLIAMS ^ in their Antelope Island mark-recapture 
experiment. 

One set of data on which the model was not based does, however, still remain. That data is 
the measured initial water reserve for a number of species. It has the makings of a test for 
consistency with reserves. In a perfect world, the measured, critical water loss contour should 
correspond to that of this model's, median or 50% emergence contour. G. morsitans is of 
obvious interest as the species on which the model is based. 




90 
70 
50 
30 
10 



Figure 8: Computed pupal emergence (left) and water loss (right) for G. morsitans. 
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Figure 9: Computed pupal emergence (top left) and water loss (top right) for G. brevipalpis; 
computed pupal emergence (bottom left) and water loss (bottom right) for Glossina pallidipes. 
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Figure 10: Computed pupal emergence (top left) and water loss (top right) for G. palpalis; 
computed pupal emergence (bottom left) and water loss (bottom right) for G. swynnertoni. 
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The sensitivity of G. brevipalpis and G. palpalis pupae to dehydration makes these species 
arguably the most challenging tests, as well as of particular interest to this work. G. brevipalpis 
is, furthermore, of national interest to those who funded this work and G. palpalis is known to 
feed on humans in West Africa (SOLANO ifTOl ). it being a major culprit in the spread of human 
trypanosomiasis. In both cases, information on the respective 4th instar excretions is available. 

The critical water reserve is also known for two other species which are expected to be less 
challenging as tests and for which G. morsitans-type, 4th instar excretions are expected to 
suffice. The water loss associated with the 4th instar excretions only becomes relevant close 
to the dew point (a part of the domain in which we have little interest) and it is, furthermore, 
thought to be too small to be of any real consequence. 

G. brevipalpis G. morsitans G. pallidipes G. palpalis G. swynnertoni 
18.7mg 8.8mg lO.Smg 7.7mg 8.5mg 

Table 3: Initial water reserves after BURSELL [[H. 

Given these initial water reserves, the position of the critical water loss contour should corre- 
spond to that of the median or 50% emergence contour for the model to work (Figures [8] to 
[TOl) . In this regard, it is worth noting that the G. mor^/ton^-based puparial duration used in the 
model could be as little as 83% of the G. brevipalpis puparial duration and the correspondence 
observed in the Figure [9] result should therefore not be as good as it is. 



9 Conclusions 

The results are certainly adequate to conclude proof-of-concept for the model (at very least). 

High transpiration rates are a consequence of high temperatures and unfavourable humidities. 
They lead to a dehydration of the tsetse pupa which can be fatal. Although the diametrical 
opposite is true of transpiration rates at low temperatures, metabolic processes are slowed, the 
puparial duration becomes too long and the cumulative effect of transpiration can be just as 
fatal. The pupal development of the tsetse fly may be divided into four successive stages: The 
third instar, the fourth instar, the sensu strictu pupal stage and the pharate adult. Water loss 
rates adhere in some semblance to this timetable. 

It would appear that the similarities between the various Glossina species are better than ex- 
pected so far as pupal water loss is concerned. The correspondence between actual measured, 
critical, water losses and water losses calculated based on Assumption [3l is profound in most 
cases. No discemable error exists for some species e.g. G. pallidipes and G. swynnertoni. It 
would therefore appear that most tsetse species behave fairly similarly in the ground, that they 
actively pursue a strategy to minimise water loss for the majority of modem habitats and have 
hydrational mechanisms preventative of desiccation. 
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G. Morsitans Daily Mortality 



G. Austeni Daily Mortality 
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G. Palpalis Daily Mortality 



G. Brevipalpis Daily Mortality 
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Figure 11: Computed daily mortality for G. morsitans (top left), G. austeni (top right), G. 
palpalis (bottom left), G. brevipalpis (bottom right). 
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Ratio of Females to Males (G. Morsitans) 
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Figure 12: Computed ratio of female to male emergent for G. morsitans (left) and G. bre- 
vipalpis (right). Note that this is a minimum scenario, owing to the lack of availability of 
any gender based data on pupal masses and surface areas (these results are based on puparial 
duration alone). 



The Glossina genus very likely does derive from a common, tropical, rain-forest dwelling 
ancestor, adjusted to moist, warm climates and dehydration is a challenge to pupae. The fact 
that loss rates in BURSELL |[1] were originally determined as fractions of initial pupal masses, 
rather than per unit of surface area, turned out to be a windfall rather than a criticism (defended 
in BuRSELL O). The necessary information to convert between the two turns out to facillitate 
the adaption of the G. morsitans model to other species. 

A strong school of thought is of the opinion that pupal and teneral mortality due to dehydration 
are either irrelevant or can be assumed constant. Do daily, pupal mortality and water loss 
justify this work? A cursory inspection of Figure [HI is suggestive of what the origins of 
such an argument might be in the case of the Morsitans group, even for G. austeni^ . Notice, 
however, that even for a hardy fly such as G. morsitans, its prospects deteriorate rapidly once 
out of favourable habitat. Even for the morsitans group, daily pupal mortality is neither linear, 
nor a function of temperature alone. When it comes to some members of ihtfusca and palpalis 



^G. austeni is thought to have secondarily invaded moist forest areas (it is small, with all the lack of hydrational 
inertia that a high surface area to volume ratio implies). 
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groups, however, that pupal mortality due to dehydration is both relevant and palpable is beyond 
contention. The effect of humidity is profound. Humidity defines habitat. This is despite the 
fact that water loss and any consequent pupal mortality are also very different things (one 
expects water loss to culminate in, and ultimately take its toll on the teneral). 

Two gender selection effects are evident based on puparial duration alone. Under adverse 
conditions more females than males emerge. A less significant, male-selective phenomenon 
occurs close to dewpoint. Both phenomena are expected to be enhanced should the heavier 
female mass and lower relative surface area be taken into account (no data could be found). 

The prognosis for the simplistic experimental model would seem to be better than expected 
(given that issues such as inferior quality pupae, differing puparial durations and the shortage 
of statistically significant data can be corrected at some stage). 

The possibility that plots presented in this work may point to the fact that breeding sites for 
some species could be predicted to be very much confined in the dry season, cannot be ruled 
out. These would be obvious places in which to concentrate traps and one immediate applica- 
tion of this work. 
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Addendum 



10.1 The Issue of Sub-Standard Laboratory Pupae 

In Burs ELL it is somewhat heuristically argued that the laboratory pupae in question were 
too small and that all emergence curves should therefore be displaced 10% to the left (the right 
in those graphs). An alternative argument can also be entertained: To bring pupae that are 

10% too small up to size would mean multiplying volume by or each dimension by ^/^. 
It would not be unreasonable to assume reserves are volume dependent whereas water loss is 
dependent on surface area. 

Suppose $ is the initial starting reserve of the pupa and that g{T, t) is some complicated loss 
rate function. If emergence is dependent on some critical water loss, 

$-g{T,t) s {100- h) = 0, (18) 

for the larger, wild puparium, one then has 

2 

^$-9(r,t) (H)%(ioo-;/) = 0, (19) 

where H is the humidity pertaining to a larger, wild puparium. If equation [18] was satisfied for 
h, then equation [19] will be satisfied for the new humidity, H, when 



Thus, replacing the argument in the old emergence relation, E(h), with 

{H + 3.57) 



is as good as the new one. 



E 



H + 3.57 



so far as puparial transpiration is concerned. Although puparial transpiration rates are of the 
order of ten times bigger than pupal rates, the sensu strictu pupal period is usually of the order 
of ten times as long. In other words, the pupal phase is just as relevant, the same reasoning as 
above can be applied to the historical dependence (the coefficients for equation [6] suggest that 
historical effects are as important as the prevailing humidity). Dependence on H is quadratic. 

Of course, observed emergence is not linearly dependent on humidity, as the relationships 
depicted in Figure [7] show. If, however, one is talking about the relevance of size to straightfor- 
ward water loss (and water loss being the predominant determinant of emergence), not variation 
within the species, one would question a simple, sideways shift of 10%. 
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10.2 Does Generalizing the Morsitans Based Model to other Species Work? 

On the face of it, Assumption [3] is certainly the most tenuous. How valid is it? Does such a 
simplistic approach work? Both the examples which follow suggest that the model's general 
puparial rate is slightly on the low side for the first few hours, something we already know to 
be true. 



Example 1 

Consider the conversion of first day G. morsitans transpirational values to G. brevipalpis. Then 
dk 

-^{h,T) = 0puparium(/i) ^puparium(T) 5pup^,um = I-IO X lO'^ X 10.3 = 11.4 X lO'^^h'^ 

This answer is, of course, in dimensionless, G. morsitans pupal masses. To convert to a fraction 
of G. brevipalpis pupal mass h~^: 

11.4 X lO-^h-i X Uhl!^^^ = 11.4 X 10-\-' X |i = 4.51 x IQ-^h"^ 

'^brevipalpis ' " 

Compare this with the measured value of 4.79 x 10"^ ± 0.05 x lO"'^ pupal masses h~^ cited in 
BURSELL This constitutes a 6% error. 



Example 2 

Consider the conversion of first day G. morsitans transpirational values to G. palpalis. Then 
dk 

-^{h,T) = (/)puparium(/i) ^puparium(T) ^p^parium = 1-10 X lO'^ X 2.54 = 2.80 X lO'^^h-^ 

This answer, again, is in dimensionless, G. morsitans pupal masses. To convert to a fraction of 
G. palpalis pupal mass h"^: 

2.80 X 10-\~' X = 2.80 X lO-^h-^ x — = 2.70 x lO-^h-^ 

'^palpalis 32 

Compare this with the measured value of 3 x 10^^ ± 0.4 x 10~^ pupal masses h^^ cited in 
BuRSELL [1|. The 10% error is within the measured limits and tolerable by 'engineering 
standards'. 



